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1 Introduction 

Our main goal in this work is to study magnetized neutron stars by using a fully 
general—relativity approach presented in the LORENE packagli]. Here we have 
adopted a non-uniform magnetic held prohle which depends on the baryon density. 
This prohle has been used in many previous works and seems to be a good choice 
to explore maximum ehects of the internal magnetic held in these objects. Equally 
important, the magnetic held treated here is poloidal and axisymmetric. The prelim¬ 
inary results show that stars endowed with a strong magnetic held will be deformed 
and the mass somewhat increased. 


2 Formalism 

The choice of the coordinates in General Relativity is crucial not only to write the 
gravitational equations in an advantageous form, but also it can make the problem 
easier to solve numerically. In the present case, due to the symmetry of the problem, a 
polar—spherical type coordinates is chosen, namely, the Maximal—Slicing—Quasi—Isotropic 
coordinates (MSQI) [l|. The metric in the MSQI coordinates is written as: 

ds^ = —N'^dt^ -\- sin^ Oidf — — ^{dr"^ -\- r'^dOY (1) 

with N'f’ir^d) the shift vector, N{r,9) the lapse function and A and B are functions 
of r and 6. Details of the gravitational equations can be found in the references [2|-[^. 

The energy momentum tensor of the system reads: 

Taf} = {e+ p)Uay + POaf} + ^ (2) 

where e is the energy density and p the pressure of the fluid. The second term is the 
electromagnetic contribution and we are not taking into account the magnetic held 

^http://www.lorene.obspm.fr 
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in the equation of state. Work along this line is in progress. Details of the equation 
of state used in this work can be found in the reference j^. 

The magnetic held measured by the Eulerian observer is given by [^: 

Ba = = (3) 

fc] 1 n\ 

\ ’ A^Br"^ sin 9 89 ’ A^B sin 9 dr ' / ’ 

where eaji'ycr is the Levi - Civita tensor related to the metric and the four 
velocity of the Eulerian observer. Assuming that the matter inside the star has 
inhnite conductivity, the electric held measured by the comoving observer must be 
zero. 

The stress—energy tensor of the magnetic held (second term in eq. |2]) is: 

IT = T (4) 


from which one can obtain the sources of the gravitational helds. The energy density 
reads: 

r.EM 1 + /c, A \2 




sin^ 6* ^ ’ 


(5) 


with U the velocity of the huid in the cp direction. The momentum density can be 
written as: 

^ ^ (dA^f (6) 


jEM _ 

^ An A^r sin 9 


and the stress 3-tensor components are given by: 


oEMr 

O yA 


1 1-f/^ 1 (dA^V 

Stt sin^ 9 \ / r'^ \ 89 ) 


(7) 


qEM6 
^ 6 


1 l-D^ 1 (dA^y 
Sir A^r"^ sin^ 9 \ 89 ) 



( 8 ) 


i^EA4(p 


1 1 - [/2 
Stt A®r2 sin^ 9 


(dA^f 


(9) 


All the sources term depend on the magnetic vector potential A^. Besides the 
above dehnitions, others important quantities measured by Eulerian obeserver are 
the circunferencial radius Rdrc 


TT TT 

Rcirc ^ (^^9 5 


( 10 ) 
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being r^q the coordinate equatorial radius, the total gravitational mass of the star 


/ NA^ 2 

+ S' + sin Odrdddcf), 

and an ellipticity which quantihes the apparent oblateness reads 


e = 



( 11 ) 


( 12 ) 


3 Magnetic field profile 

Magnetars are know to possess strong magnetic helds which can be estimated from 
the the period and the period derivative of the star. These highly magnetized neutron 
stars are believed to have surface magnetic held of order of 10^^ — 10^^ G . Besides, 
the virial theorem states that the magnetic held can be still much higher inside the 
star than at the surface. Following the references [8 — 16], a non-uniform magnetic 
held which can depend on the density is parametrized as: 


=B, + Bo 

1 — e y^oJ 




The parameters are the nuclear density at saturation no, the barion denstity of matter 
nfc, the magnetic held on the surface Bg, and a parameter that controls the magnetic 
held at the center Bq. We set the parameters (3 and 7 to j3 = 0.01 and 7 = 2.0. Other 
choices for these parameters are possible, but qualitatively the main conclusion remain 
the same. Usually the value of the magnetic held at the center is about 70% of So- 
Based on the equations dehned in the previous sections, in order to construct our 
models we need to have an expression for the vector magnetic potential A^. As we 
can note in flT^ there is no way to get the magnetic potential directly, since the 
magnetic held is dehned as a function of the baryon density ni,. Alternatively, if we 
suppose a constant magnetic held in the z direction B = Bz, the magnetic vector 
potential is given by: 

7^ yc B 1 

A^ = — - —Ri-rSsin00. (14) 

Our hrst approximation is to say that the magnetic held B in equation flTT|) is 
locally constant and described by the equation flT^ from which we can easily see that 
the magnetic held is conhned inside the star, reaches its maximum value at the center 
of the star. 

To illustrate the approach we show on the Figure 1 isocontours of the magnetic 
held strength, which linies lie on the surfaces A^ = const. For a specihc choice of 


3 







Magnetic field: 



X [km] 

Figure 1: Isocontours of the magnetic field strength in the (x, y) plane 

central energy density we have obtained a mass of M = 1.31 Msun and a circunference 
radius of Rdrc = 11-34 km, whose result are very close to the spherical case and 
without magnetic held. The ellipticity was found to e = 0.032. This tiny effect is due 
to low value of the magnetic held, namely, B = 3.5 x 10^^ G at the center, which is not 
enough to deform the star. In order to investigate how the mass and the ellipticity, 
and therefore, the deformation, change with the magnetic held we have evaluated 
models for diherent values of Bq. 

Figure 2 shows the change in the mass, while the Figure 3 presents the ellipticity 
as a function of the magnetic held. 

From these hgures, both the mass and the ellipticity increase with the magnetic 
held whose maximum value corresponds to a value of 1.64 x 10^® G at the center of 
the star. The corresponding ellipticity is found to be e = 0.146, showing that the 
deformation plays an important role in the global properties of the star. According 
to the Fig.2 the increasing in the mass is about 7%, less than previous calculation 
by solvin g T OV equation, and higher than the 2% predicted by the perturbative 
approach |l6| . 
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Figure 2: Curve of the gravitational mass M as a function of central magnetic field. 
0.16 

0.14 

0.12 

0.1 

I . 0.08 
■© 

0.06 

0.04 

0.02 

0 

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

Bc(lo'^G) 

Figure 3: Curve of the ellipticity e as a function of the central magnetic field. 

4 Conclusion 

We have used the Lorene package for computing perfect fluid magnetized stars in 
general relativity with the inclusion of a non-uniform magnetic held prohle. We have 
then used the code to see the effect of the magnetic held strength on the deformation 
and also on the mass of the star. The maximum magnetic held found at center is 
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around 1.64 x 10^®G for a fixed stellar mass of M = 1.40 The results were 

obtained for the APR equation of state. A study with a more sophisticated equation 
of state including magnetic effects and more realistic held conhgurations is in progress. 
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